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Abstract 

Different models of random walks on the dual graphs of compact ur- 
ban structures are considered. Analysis of access times between streets 
helps to detect the city modularity. The statistical mechanics approach 
to the ensembles of lazy random walkers is developed. The complexity of 
city modularity can be measured by an information-like parameter which 
plays the role of an individual fingerprint of Genius loci. Global struc- 
tural properties of a city can be characterized by the thermodynamical 
parameters calculated in the random walks problem. 



PACS codes: 89.75.Fb, 89.75.-k, 89.90.+n 
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1 Introduction to city networks studies 

Studies of urban networks have a long history. Many researches have been 
devoted to the optimizations of transportation routes and power grids, to the 
predictions of traffic flows between the highly populated city districts, to the 



*The Alexander von Humboldt Research Fellow at the BiBoS Research Center 



1 



investigations of habits and artefact exchanges between distanced settlements in 
historical eras. In the most of them, relations between different components of 
urban structure are often measured along streets, the routes between junctions 
which form the nodes of an equivalent planar graph. The graph-theoretic princi- 
ples had been applied in [1] in order to measure the hierarchy in regional central 
place systems, in [2] to the measurement of transportation networks. The use 
of graph-theoretic view and network analysis of spatial systems in geographic 
science had been reviewed in It is interesting to mention that graphs have 
been widely used to represent the connectivity between offices in buildings |4J 
and to classify various building types in [5]. 

In all these studies, the traffic end points and junctions had been treated 
as nodes, and the routes had been considered as edges of some planar graphs. 
Being embedded onto the geographical and economical landscapes, these pla- 
nar graphs bare their multiple fingerprints. Among the main factors featuring 
them are the high costs for the maintenance of long-range connections and the 
scarce availability of physical space. Spatial networks differ from other complex 
networks and call for the alternative approaches to investigate them [6j. 

While studying the motifs and cycles in the complex networks, a compara- 
tive analysis of the original graph and of its randomized version is used. If a 
motif is statistically significant, it appears in the real network much frequently 
than in the randomized versions of the graph 7 J. However, in the planar city 
street patterns, its randomized version is not of significance since, first, it is 
surely a non-planar graph due to the randomness of edge crossings and, second, 
the long-range connections which inevitably present in the random graphs in 
abundance are extremely costly in the real cities [BJ. In [5], it has been pro- 
posed to compare the city street patterns with the grid-like structures that is 
indeed useful essentially for the regular urban development. 

It had been formulated in the classical essay of T. Harold Hughes [9] that 
the accessibility of important city objects for the vehicle traffic and pedestrians 
is always the chief factor in regulating the growth and expansion of the city. A 
broad, simple scheme of main traffic lines gives a sense of connectedness and 
unity to the various parts of the city and links up country and town. In [10J [11] , 
a significant correlation between the topological accessibility of streets and their 
popularity, micro-criminality, micro-economic vitality and social liveability had 
been established. In the traditional representation of space syntax based on 
relations between streets through their junctions, the accessibility or distance is 
associated with points or junctions. 

We would like to mention that the issues of global connectivity of finite 
graphs and accessibility of their nodes are the classical fields of researches in 
graph theory. They are studied by means of certain dynamical processes defined 
on the graphs. In particular, in order to reach "obscure" parts of large sets and 
estimate the probable access times to them, random walks are often used [12]. 
There are a number of other processes that can be defined on a graph describing 
various types of diffusion of a large number of random walkers moving on the 
network at discrete time steps [14]. In all such processes we have deal with 
discrete time Markov chains studied in probability theory. Markov chains have 
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the property that their time evolution behavior depends only upon their current 
state and the state transition properties of the model. 

At each time step every walker moves from its current node to one of the 
neighboring nodes along a randomly selected link. The metric distance between 
the nodes are of no matter for such a discrete time diffusion process and then the 
focus of study is naturally shifted from the original problem of traditional space 
syntax to the dual one based on relations between the streets which themselves 
are treated as nodes. The distance between two streets, in such a representation, 
is a distance in the graph-theoretic sense. The dual graphs of a geographic net- 
work comparable in its structure with other complex networks are irrelevant to 
neither distance, nor the physical space constraints. The relations between the 
traditional space syntax representation based on the relations between streets 
through their junctions and the dual representation that is a morphological rep- 
resentation of relations between junctions through their streets was studied in 
details in [13]. The dual city graphs had been developed and studied within 
the concept of space syntax |10j . The key characteristics in space syntax is that 
precedence is given to linear features such as streets in contrast to fixed points 
which approximate locations [11]. The dual city graphs have been discussed 
recently in [15] . In [16j . they have been called the information city network. 

In [17], while identifying a street over a plurality of routes on the city map, 
the "named-street" approach had been used, in which two different arcs of the 
original street network are assigned to the same street ID provided they have 
the same street name. The main problem of the approach is that the meaning of 
a street name could vary from one district or quarter to another even within the 
same city. For instance, the streets in Manhattan do not meet, in general, the 
continuity principle rather playing the role of local geographical coordinates. 

Being interested in the statistics of random walks, we generalized the ap- 
proach used in [17] in a way to account the possible discontinuities of streets. 
Namely, we assign an individual street ID code to each continuous part of it 
even if all of them share the same street name. The dual graph is constructed 
by mapping edges coded with the same street ID into nodes of the dual graph, 
and intersections among each pair of edges in the original graph - into edges 
connecting the corresponding nodes of the dual graph like it has been done in 

HZ!. 

In [15] , the Intersection Continuity principle (ICN) different from our identi- 
fication approach has been used: two edges forming the largest convex angle in a 
crossroad on the city map are assigned the highest continuity and therefore are 
coupled together acquiring the same street ID. The main problem with the ICN 
principle is that the streets crossing under the convex angles would artificially 
exchange their identifiers that is not crucial for the study of degree statistics, 
but makes it difficult to interpret the results on random walks and detect the 
dynamical modularity of the city. 

It is also important to mention that the number of street IDs identified 
within the ICN principle usually exceeds substantially the actual number of 
street names in a city. In [6] [HI [25l [26] [27] , the degree statistics and various 
ccntrality measures for the data set of the square mile samples of different world 
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cities had been investigated. However, the decision on which a square mile would 
provide an adequate representation of a city is always questionable. 

In this paper, we use an alternative strategy investigating the spectral prop- 
erties of random walks defined on the dual graphs of compact city patterns 
bounded by the natural geographical limitations. The reason we consider the 
compact urban domains is twofold. First, it allows us to avoid the problem of a 
" square mile" and, second, the compact urban domains had been usually devel- 
oped "at once", in accordance to certain architectural principles, their partial 
redevelopment had been occasional and rear, so that they can be considered 
typical. 

We have studied the Markov chains defined on the six undirected dual graphs 
corresponding to the different street and canal urban structures. Two of them 
are situated on islands: Manhattan (with an almost regular greed-like city plan) 
and the network of Venice canals (imprinting the joined effect of natural, po- 
litical, and economical factors acted on the network during many centuries). 
We have also considered two cities founded shortly after Crusades and devel- 
oped within the medieval fortresses: Rothenburg ob der Tauber (the medieval 
Bavarian city preserving its original structure from the XHI-th century) and 
the Bielefeld downtown (Altstadt Bielefeld) composed of two different parts: 
the old one founded in the XHI-th century and the modern part subjected to 
the partial urban redevelopment at the end of XlX-th century. To supplement 
the study, we have investigated the canal network of the city of Amsterdam. Al- 
though it is not actually isolated from the national canal network, it is binding 
to the delta of Amstel river forming a dense canal web showing a high degree 
of radial symmetry. The scarce of physical space is among the most important 
factors determining the structure of compact urban patterns. The general in- 
formation on the dual graphs of compact urban street and canal patterns that 
we have studied is given in Tab. 1. The spectral properties of finite Markov 
chains defined on these dual graphs are compared with those of a model exam- 
ple: a hypothetical village extended along one principal street and composed of 
N blind passes branching off it. 

A square mile of the New York city grid and a square mile pattern of Venice 
streets array have been discussed recently in [HI [25l ESI [27], however, up to 
our knowledge, the canal patterns have never been subjected to the network 
analysis. The navigation efficiency in the Manhattan streets has been studied 
in [TB]. 

It is worth to mention the importance of implemented street identification 
principle for the conclusion on the degree statistics of dual city graphs. The 
comparative investigations of different street patterns performed in [6] [15] im- 
plementing the ICN principle reveal scale-free degree distributions for the ver- 
tices of dual graphs. However, in [T7] it had been reported that under the 
street-name approach the dual graphs exhibit the small-world character, but 
scale-free degree statistics can hardly be recognized. The results on the proba- 
bility degree statistics for the dual graphs of compact urban patterns analyzed 
in accordance to the street identification principle that we have described above 
are compatible with that of [T7] . Compact city patterns do not provide us with 
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sufficient data to conclude on the universality of degree statistics. It is remark- 
able that the probability degree distributions for the dual graphs correspondent 
to the compact city patterns are broad and have a clearly expressed maximum 
and a long right tail. The presence of a noticeable maximum in the probability 
degree distributions indicates that the structures of compact urban patterns are 
usually close to a regular one and that there is the most probable number of 
junctions an average street has in a given city. The long right tails of distri- 
butions correspond to the highly connected nodes of dual graphs, just a few 
"broadways", embankments, and belt roads crossing many more streets than 
an average street in the city. To give an example, in Fig. 1, we have displayed 
the log-log plot of fractions of streets n(s) via the number of their junctions s 
in Manhattan. These number are shown by points and the solid line is for the 
relevant cumulative distribution P c (s) — J2^'=s n ( s ') 0- 

2 Random walks on the dual city graphs: the 
description of models and the sketch of results 

We consider a connected graph G = (V, E) with \V\ — N nodes and \E\ = m 
undirected edges specified by its adjacency matrix A such that Ajj = 1 if the 
node i is connected to j and Aij = otherwise. A random walk starting at a 
node vo G V and traversing a sequence of random nodes {vt} as t = 0, 1, ... is a 
Markov chain characterized by the matrix of transition probabilities T = DA, 
in which D is the diagonal matrix of inverse vertex degrees, Dij = k^ 1 5ij, where 
ki = deg(j), the degree of vertex i in graph G (the number of junctions a street 
shares on the city plan). The transition matrix T meets the normalization 
condition, Y^j^ij = 1- We denote by 7r* € M. N the distribution of vt in the 
Markov chain, 7r* = Pr(v t — i). Then the rule of the walk can be expressed by 
the simple equation 

TT^ 1 = TV (1) 

where TT is the transposed transition matrix. In the present paper, we discuss 
only the undirected dual city graphs with symmetric adjacency matrices, Aij — 
Aji. Nowadays, certain driving directions are specified for the streets by the 
traffic regulation polices, so that the relevant transportation lines appears to 
be directed. We do not consider them, limiting our present study only to the 
network of pedestrian access. The distribution of a current node in the random 
walk defined by (JTJ) on the undirected graphs after t steps tends to a well defined 
stationary distribution Tii = fc,/2m (which is uniform if if the graph is regular) 
that is a left eigenvector of the transition matrix T, belonging to the largest 
eigenvalue 1 [32], [IB]- The spectrum of problem {T]) is contained in the interval 
[-1,1]- 

In Sec. 3, we calculate and analyze the expected number of steps a random 
walker starting from node i makes before node j is visited (the access time) for 
all pairs of streets (canals) in the compact urban patterns. The properties of 
access times can be used in order to estimate the accessibility of certain streets 
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and districts by random walkers starting from the rest of city. In particular, 
the access times allows one to introduce the equivalence classes of nodes and to 
obtain a well-defined ordering of these equivalence classes, independent of any 
reference node. The nodes in the lowest class are difficult to reach but easy to 
get out of, while the nodes in the highest class are easy to reach but difficult 
to leave. We also discuss the random target access times and the distributions 
of mean access times in the compact cities. The latter characteristics can be 
used in order to detect ghettos (the groups of dynamically isolated nodes) and 
estimate the accessibility of certain districts from the streets located in other 
parts of the city. 

The problem of random walks (Q]) defined on finite graphs can be related to 
a diffusion process which describes the dynamics of a large number of random 
walkers. The associated Laplace operator is more convenient, since its spectrum 
is positive. Indeed, the eigenvalues in fTj can be negative, so that the spectral 
moments J^i Af™ cou ld oscillate strongly for large networks. In contrast to it, 
the eigenvalues of Laplacian are positive that allows one to study its spectral 
properties by the powerful methods of statistical mechanics, in which various 
functions /(Aj) defined on the spectrum {Xi} are considered. 

The diffusion process is defined by the expectation number of random walkers 
n £ M. N and described by the equation 

n = Ln, (2) 

in which the scaled Laplacian L (symmetric) is defined by 

L=l — D -U2 TD l/2 (3) 

where 1 is the unit matrix. Let us note that L is related to the transition 
matrix T in Q by T = D 1 / 2 (l - L)_D -1 / 2 . All eigenvalues of L belong to the 
interval A Q £ [0, 2], the eigenvectors are normalized, |n( Q )| = 1, and orthogonal, 
(n( Q )n(«) = 5 a(3 . 

The linear equation @ supplied with an initial condition n is solved by the 
function n* = Qexp(— tJ)(Q^ 1 n ), in which J is a block diagonal matrix (the 
Jordan canonical form of L) , Q is the the transformation matrix corresponding 
to the Jordan form. Given the eigenvalue A Q with the multiplicity m a > I, the 
relevant contribution to n' is given by 

m a /m a -l ^\ 

exp(-a Q )^ ^ c <+i7| ) u fe ( 4 ) 

k=Q \ 1=0 V 

where is the fc-th vector of the ma-dimensional eigenspace of degenerate 
eigenvalue X a and c/ is the l—th component of the transformed vector of initial 
conditions, c = Q^no. It follows from ([4]) that in the presence of degenerate 
modes with m a ;§> I the relaxation process can be delayed significantly. More- 
over, if X a < I the expectation numbers of walkers in a certain city modules 
would substantially increase in a lapse of time. 
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The analysis of spectral properties of Laplacian operator allows for detec- 
tion of fine-scale dynamical modularity of city which cannot be seen from the 
transition matrix T and provide us a tool for the estimation of entire city sta- 
bility with respect to occasional cuts of certain transportation lines breaking it 
up into dynamically isolated components. In Sec. 4, we consider the modes of 
diffusion process defined on the dual city graphs and compute the coefficients 
of linear correlations between densities of random walkers that flow along the 
edges of dual city graphs in the discrete steps. Specifying the certain correlation 
threshold, we detect the groups of nodes traversed by the essentially correlated 
flows of random walkers. The structures of essentially correlated flows and their 
appearance are the individual characteristics of a city. We measure the quantity 
and extension of clusters of nodes sharing the essentially correlated flows of ran- 
dom walkers by an information parameter reflecting the complexity of random 
walk traffic. 

In Sec. 5, we discuss the statistical mechanics of so called lazy random walks 
specified by the parameter < p < 1. In the model of lazy random walks, 
an agent located at a node v moves to a neighboring node with probability 
Pky 1 , but rests modeless with probability (1 — /?). We compute the well-known 
thermodynamical quantities (the internal energy, entropy, the free energy, and 
pressure) describing the macroscopical states of an ensemble of " lazy" random 
walkers on the dual graphs of compact urban patterns. We provided the detailed 
interpretation for each thermodynamical parameter in the context of random 
walks. 



3 Access times in the compact city structures 

The important characteristics of random walks defined on finite graphs is the 
access time H^, the expected number of steps before node j is visited, starting 
from node i |12j . The elements of matrix are computed for each pair of 
nodes i,j following the formula: 

N / 2 \ 

in which /ii = l>/i2>---MJV> — 1 are the eigenvalues and ip s are the relevant 
eigenvectors of the symmetrized transition matrix D~ 1 / 2 TD 1 / 2 . 

The access time from i to j may be different from the access time from j to 
i, Hij ^ Hji, even in a regular graph. A deeper symmetry property of access 
times for undirected graphs was discovered in [19] . 

+ Hjk + Hki = Hik + Hkj + Hji (6) 

for every three nodes in G. This property allows for the ordering of nodes in 
the graph with respect to their accessibility for the random walkers. It has 
been pointed out in [12] that the nodes of any graph can be ordered so that if i 
precedes j then Hij < Hji. This ordering is not unique, but one can partition 
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the nodes by putting i and j in the same equivalence class if H^ — Hji and 
obtain a well-defined ordering of the equivalence classes, independent of any 
reference node. The nodes in the lowest class are difficult to reach but easy to 
get out of, while the nodes in the highest class are easy to reach but difficult to 
leave. If a graph has a vertex-transitive automorphism group, then = Hji 
for all i,j £ G. The random target identity [18], 

^irjHij = Const, (7) 

3 

states that the expected number of steps (the random target access time, r) 
required to reach a node randomly chosen from the stationary distribution -k is 
a constant, independent of the starting point of the given graph G. The values 
of random target access times grow with the size of graphs and are very sensitive 
to their structures. Their values for the compact urban structures are given in 
the last column of Tab. 1. 

The properties of access times can be used in order to estimate the accessi- 
bility of certain streets and districts by random walkers starting from the rest 
of city. Computations of access times to the streets in the studied compact 
urban structures convinced us that for any given node i the access times to it, 
Hij, change with j inferentially in comparison with their typical values, and 
therefore, the mean access time, 

1 ™ 

i=i 

can be considered as a good parameter for estimating the accessibility of a street 
by random walkers. Distributions of mean access times, a/,, in the city can be 
consider as estimations of its connectedness. In particular, it helps to detect the 
ghettos, the groups of streets almost isolated (in the dynamical sense) from the 
rest of town. In Fig. 2, we have presented together the distributions of mean 
access times to the streets in Manhattan (dashed line) and Rothenburg o.d.T. 
(solid line). The distribution for Rothenburg o.d.T. exhibits a local maximum 
at relatively long access times (rj 300 steps) indicating that there is a number 
of low accessible streets in the town. The distributions of mean access times to 
the canals in Amsterdam (dashed line) and Venice (solid line) is presented on 
Fig 3. 

The distribution of mean access times in the downtown of Bielefeld is of es- 
sential interest since it comprises of two structurally different parts (see Fig. 4. a). 
The part "A" keeps its original structure (founded in XIII-XIV cs.), while the 
part "B" had been subjected to the partial redevelopment in the XlX-th century 
(Fig. 4. a). It is important to mention that the city districts constructed in accor- 
dance to different development principles and in different historical epochs can 
be easily visualized on the dual graph of the city. In Fig. 4.b, we have shown the 
3D-representation of the dual graph of the Bielefeld downtown. The (xi, Hi, Zi) 
coordinates of the i— th vertex of the dual graph G in three dimensional space are 
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given by the relevant i— th components of three eigenvectors tS 2 \ u^^and 
of the adjacency matrix Ac (which does not coincide with the transition matrix 
T). These eigenvectors correspond to the second, third, and fourth largest (in 
absolute value) eigenvalues of Ac [20] . The 3D-dual graph of Bielefeld displays 
clearly a structural difference between "A" and "B" parts: in 3D- represen- 
tation, the relevant subgraphs are located in the orthogonal planes (Fig. 4.b). 
Sometimes other symmetries of dual graphs can be discovered visually by using 
other triples of eigenvectors if the number of nodes in the graph is not too large. 

In Fig. 5, we have displayed the distributions of mean access times h to the 
streets located in the medieval part " A" starting from those located in the same 
part of Bielefeld downtown, from "A" to "A" (solid line). It has been computed 
by averaging over i, j 6 A in ©. The dashed line presents the distribution 
of mean access times to the streets located in the modernized part " B" starting 
from the medieval part "A" (from "A" to "B", i G A and j £ B). One can see 
that in average in takes longer time to reach the streets located in " B" starting 
from " A" . The similar behavior is demonstrated by the random walkers starting 
from "B" (see Fig. 6): in average, it requires longer time to leave a district for 
another one. Study of random walks defined on the dual graphs helps to detect 
the quasi-isolated districts of the city. 

4 Dynamical modularity in the compact city 
structures 

Dynamical modularity is a particular division of the graph into groups of 
nodes on which the certain modes of diffusion process are localized. It can be 
detected by analyzing spectral properties of the relevant dual graph that is of the 
spectrum of some differential operator defined on it |21j . In particular, for the 
undirected graphs, it is convenient to consider differential elliptic self-adjoint 
operators with a positive spectrum. We have studied the Laplace operator 
defined by ©. 

The eigenvalues of Laplace operator ([3]) along with the continuous approx- 
imations of their densities for the dual graphs of compact urban patterns are 
shown in Figs. 7-9. The similarity of spectra allows us to divide the studied 
compact urban patterns into three categories: i) the medieval cities: Bielefeld 
(Fig. 7) and Rothcnburg o.d.T.; ii) the canal patterns: Venice (Fig. 8) and 
Amsterdam; hi) the spectra with a highly degenerated eigenmode (A = I): 
Manhattan (Fig. 9). 

It is important to note that the densities of eigenvalues for the compact 
urban patterns differ dramatically from those computed for the classical random 
graphs of Erdos-Renyi model [H] [55], deviate form the semicircular law, and 
the densities found for the scale-free random tree- like graphs in [24] . In [28] , the 
density of eigenvalues for the Internet graph on the Autonomous Systems (AS) 
level had been computed. The Internet spectrum on the AS level is broadly 
distributed with two symmetric maxima and similar to the eigenvalue density 
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of random scale-free networks. In contrast to them, the spectral density for 
the compact city samples are bell shaped and tend to turn into a sharp peak 
localized at A = 1 due to the highly degenerate 1-mode which would score a 
valuable fraction of eigenvalues (48% for Manhattan). 1-modes come in part 
from pairs of streets of minimal connectivities branching of either "broadways" 
or belt roads. These structures are overrepresented in the compact city patterns. 

The slowest modes of diffusion process ((2]) allow one to detect the city mod- 
ules characterized by the individual accessibility properties. Indeed, these ran- 
dom walks do not describe the properties of an actual city traffic, since the 
approach does not concern the lengthes of streets referring uniquely to the topo- 
logical properties. However, the modularity plays an important role in any local 
dynamical process taking place in the city including its traffic conditions. The 
primary feature of the diffusion process in the compact urban patterns is the flow 
between the dominant pair of city modules: the "broadways" and the relatively 
isolated streets remote from the primary roads. 

Due to the proper normalization, the components of the eigenvectors n( a ) 
play the role of the Participation Ratios (PR) which quantify the effective num- 
bers of nodes participating in a given eigenvector with a significant weight. 
This characteristic has been used in [28] and by other authors to describe the 
modularity of complex networks. However, under the abundance of a highly de- 
generate mode PR is not a well defined quantity (since the different vectors in 
the eigensubspace corresponding to the degenerate mode would obviously have 
different PR). 

As time advances the distribution of random walkers approaches a steady 
state n°° oc fcj, in which all diffusion currents are balanced. It corresponds 

to the principal eigenvector related to the smallest eigenvalue of L. The 
relaxation processes toward the steady state are described by the remaining 
eigenvectors n( a \ a > 1, with the characteristic decay times r^ a \ such that 
exp(— 1/t^ 01 ') = A Q . The second smallest eigenvalue of the scaled Laplacian ([3]) 
is related to the graph diameter, diam(G) < —In (AT — l)/ln(A2), the maximum 
distance between any two vertices in the graph. It is also related to the Fiedler 
vector [2jjj describing the algebraic connectivity of the graph. Namely, let us 
consider the components of the eigenvector n*- 2 ) corresponding to the second 
smallest eigenvalue A2 for the scaled Laplacian ([3]) defined on a connected graph 
G(V,E). Define Vi = {v G V : nf ] < 0} and V 2 = {v e V : > 0}, then the 
subgraphs induced by V\ and V% are connected. 

In general, each nodal domain on which the components of the eigenvector 
n( Q ) does not change sign refers to a coherent flow (characterized by its decay 
time 7"'°)) of random walkers toward the domain of alternative sign. The nodal 
domains participate in the different eigenmodes as one degree of freedom, and 
therefore their total number is important for detecting the dynamical modular- 
ity of city networks. It is known from J3UI that the eigenvector n'") can have at 
most a + m a — 1 strong nodal domains (the maximal connected induced sub- 
graphs, on which the components of eigenvectors have a definite sign) where m a 
is the multiplicity of the eigenvalue X a , but not less than 2 strong nodal domains 
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(a > 1) [31 j . However, the actual number of nodal domains can be much smaller 
than the bound obtained in [3T]. In the case of degenerate eigenvalues, the sit- 
uation becomes even more difficult because this number may vary considerably 
depending upon which vector from the m Q -dimensional eigenspace of degener- 
ate eigenvalue A Q is chosen. A fragment of nodal matrix for the Chelsea village 
(Fig. 10. a) in the Manhattan island is shown on Fig. 10. b: the components of 
all eigenvectors localized on the nodal domains displayed in white (black) have 
always the positive (negative) sign. 

We investigate the statistics of components of eigenvectors localized on 
a given street, {n\ } a . To uncover a fine modular structure of compact city 
patterns, we have computed the linear correlation coefficients between the lists 
of eigenvector components for all pairs of streets in a city, 



The linear correlation measures how well a linear function explains the rela- 
tionship between two data sets. The correlation is positive if an increase in the 
eigenvectors components localized on a given street corresponds to an increase in 
those related to other street, and negative when an increase in one corresponds 
to a decrease in the other. 

The study of correlations between the components of eigenvectors allows us 
for a precise recovering of all dynamical modules in a city, street by street. By 
tuning the sensitivity threshold k > 0, one can detect the groups of streets char- 
acterized by significant pairwise correlations Cy > k. In case of a degenerate 
eigenmode, the whole bunch of streets of minimal accessibility joins a correlation 
cluster at once. Investigating the size of clusters characterized by the significant 
pairwise correlations between the eigenvectors components, we have found that 
it is very sensitive to the threshold value k. Pairs of correlated streets can be 
detected for k < n c , below the critical value k c individual for each city. For 
instance, the strongest correlations in the Manhattan island (k > 0.16) are ob- 
served between the Bowery st. (Chinatown) and the 6-th Avenue (Greenwich 
Village), and between the Lafayette st. and Mott st. (Chinatown). By reducing 
k just by 1%, we immediately get many new correlated pairs and [South FDR 
Dr., Allen st.], [Greene st.(Soho), Ave. D], [Crosby st. (Little Italy), Ave. D], 
[Centre st, Elizabeth st.], [Allen st., 6-th Ave.], [Lafayette st., Ave. B], [Broad- 
way, West st.] among the others. In most cases the streets in correlated pairs 
have the same driving directions although the driving direction data has not 
been initially used in the adjacency matrices. Reducing the value of sensitivity 
threshold, one can detect the triples of correlated streets and more structurally 
complicated dynamical modules. Then the correlated clusters merge exhibiting 
a sharp phase transition into a giant connected component for essentially small 

K. 

The appearance of correlated fc-tuples and their total numbers at given k 
are by no means random and encode the important information on the city 
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connectedness. The complexity of dynamical modularity can be measured by a 
quantity of information encoded by the number of various /c-tuples of essentially 
correlated nodes at the different values of correlation threshold n in the following 
way. Since the probability that a connected correlated fc-tuple appears in a 
random labelled graph of N nodes is 



its appearance corresponds to some quantity of information I = —pk(N) log 2 pk(N). 
For a mesh of values k > 0, we had counted the total numbers of all correlated 
fe-tuples appeared in each city, A4(k), and then found the total amounts of 
information T(k) encoded by the dynamical modularity in each studied city 
pattern, 



It is worth to mention that in the absence of correlations as well as in the case 
when the all fluxes of random walkers though nodes of the dual city graph are 
correlated, no information can be encoded, and therefore 1=0. The results on 
the comparative information analysis of dynamical modularity of the compact 
urban structures are displayed in Figs. 11-13. 

In general, the quantity of information encoded by the dynamical modularity 
of ancient cities exceeds that in the modern or redeveloped ones since a number 
of correlated clusters of diverse sizes could appear in the graphs with a less 
regular structure. Information (via the threshold of correlations, n) calculated 
for Rothcnburg and Bielefeld (subjected to a partial redevelopment) are almost 
equal (Fig. 11). Both medieval cities contain a number of streets characterized 
by the relatively highly correlated traffic (they are the belt roads encircling the 
cities along their fortress walls), while the traffic along the streets close to the 
city centers appears to be less correlated, although all streets of low accessibility 
join the correlation cluster at once at some level k, then the information on the 
dynamical modularity turns to zero. 

Information profiles obtained for the Venice and Amsterdam canal networks 
look similar (see Fig. 12). However, it seems that information in a message on 
possible correlations of gondoliers traffic in Venice could be much more valuable 
than a cruise schedule along the Amstel embankments. The triggering between 
different information states displayed on Fig. 12 corresponds to the merging of 
diverse correlated clusters into the bigger correlated modules. The information 
peaks located close to k = indicate a phase transition to a giant correlated 
component which cover most of the nodes in the networks. Nevertheless, not 
all canals join the giant correlated module (since the amount of information is 
not zero as k — * 0, but, on the contrary, turns to be high). It happens because 
of the presence of an anii-correlated module in the canal network, in which the 
components of Laplacian eigenvectors are negatively correlated, dj < 0. 

The scale of the graph shown in Fig. 13 is incompatible with those of Fig. 11 
and Fig. 12 since the dynamical modularity in the Manhattan island (one of 





(10) 



k— tuples 
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the most regular city grids in the world) contains essentially less information (it 
comes primarily from the correlated clusters risen in the region of Central park) 
than that in the ancient cities. 



5 Statistical mechanics of lazy random walks in 
the compact city structures 

A number of different " measures" quantifying the various properties of com- 
plex networks has been proposed in so far in a wide range of studies in order to 
distinguish the groups of nodes and enlighten the relations between them. Some 
measures can be computed directly from the graph adjacency matrix: likelihood 
[36] . assortativity [37], clustering [38], degree centrality [33], [39], [40], betweeness 
centrality [JD], link value [UJ, structural similarity [35], distance (counting the 
number of paths between vertices) [42]. Other measures (concerned with the 
networks embedded into Euclidean space) involve the lengths of links or the true 
Euclidean distances between nodes: closeness centrality [53], [33], straightness 
centrality [44] . [45] , expansion [41] . information centrality and graph efficiency 
[45] . [46] . A good summary on the several centrality measures can be found in 
[27] and in [47] for the Internet related measures. The list of available measures 
is still far from being complete, the new measures appearing together with any 
forthcoming network model. It is also worth to add some spectral measures 
(concerned the eigenvalues of graph adjacency matrix): subgraph centralization 
[48] . subgraph centrality [34] , network bipartivity [49] and many others. 

In the present section, we study the statistics of flows of lazy random walk- 
ers roaming in the compact city patterns. The random walkers have no mass 
and do not interact with each other, so that they do not neither to contribute 
to the energy nor to the momentum transfer. Nevertheless, their flows have 
the nontrivial thermodynamical properties induced by the complex topology of 
streets and canals they flow along. The obvious advantage of statistical me- 
chanics is that the mathematical objects we introduce and all relations between 
various statistical quantities are well known in the framework of thermodynamic 
formalism. 

In the following, we use the inverse temperature parameter > which can 
be considered cither as an effective time scale in the problem (the number of 
streets a walker passes in one time step) or as the laziness parameter defined in 
Sec. 2. Since the eigenvalues A a of scaled Laplacian operator ([2]) are positive 
and bounded, one can define for them three well known spectral functions: 

1. the heat kernel (the partition function), 



N 




(11) 



converging as j3 > 0, for N — > oo 
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2. the spectral zeta function (the spectral moments) 

N 

c(s) = 5>« s , ( 12 ) 

a=l 

3. and the spectral density, 

N 

p(X) = J2S(X-X a ). (13) 

a=l 

These spectral functions are related to each other by the Laplace transformation, 

poo 

K{p) = / dX e^ A p(A), (14) 
Jo 

and by the Mellin transformation (up to the T- function) [50] . 

^ = tTT / (15) 

r(s) 7 

Furthermore, 

/oo 
P (A)A- s dA. (16) 
-oo 

Let us note that since < A Q < 2, we immediately obtain bounds for the 
spectral moments, 

C(-n) < N ■ 2 n , C(n)>N-2- n . 
The partition function K(0) meets the conditions of Bernstein's theorem [32] : 

(-l) n KW(p) > 0,/3 > 0,7i > 1, 

and therefore defines a unique non-negative measure d/i(/3) = K(j3)dj3 on R + 
providing the solutions for the Hamburger and Stiltjes moments problems, [50] . 
Concerning the list of "measures" given above, it is worth to note that, in 
principle, any of spectral functions (|11|13|) at a given temperature can be 
used as a measure of some quantities relevant to a certain lazy random walk 
model. 

The ensemble of random walkers on the graph G(V, E) can be characterized 
by the following macroscopical quantities: 1) the internal (averaged) energy, 
E = -dp\nK{(3); 2) the entropy, S = \nK{(3) + f3E; 3) the free energy, F = 
(3^ 1 InK(fl), and 4) the pressure, P = (3~ 1 d\y\ In K(/3). In the last equality, 
we have considered the differential with respect to the graph size as d\V\ = 
p(X)dX. The thermodynamics of compact urban city patterns is presented on 
the Figs. 14-22 and in Tab. 2. The collected data give an insight into the various 
aspects of averaged streets (canals) accessibility in a given city. 
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Due to the complicated topology of streets and canals, the flows of random 
walkers exhibit the spectral properties similar to that of a thcrmodynamical 
system characterized by a nontrivial internal energy (Figs. 14-16). It grows 
with temperature (the inverse parameter of lazy random walks) 0~ x (albeit still 
bounded in the interval [0, 2]). As usual, the absolute value of internal energy 
relevant to a given city cannot be precisely measured, but we can measure its 
difference in any temperature interval. In principle, the slopes of internal energy 
curves are steeper (the internal energy grows faster with f3~ 1 ) in the modern 
cities with a quite regular grid-like structures, in which any street has a relatively 
good accessibility. 

In thermodynamics, entropy is an extensive state function that accounts for 
the effects of irreversibility in thermodynamic systems. It describes the number 
of the possible microscopic configurations of the system. In the problem of finite 
random walks, its value quantifies the diversity of flows which can be detected 
in the city at a given temperature and grows with temperature. In Figs. 17-19, 
we have displayed the entropy curves via /3 _1 . It is important to mention that 
the improve of street accessibility causes a decrease of entropy growth rates for 
large temperatures (small 0). The entropy of random walker flows as well as its 
growth rate in Manhattan is ever less than in any other compact city structure 
we studied. 

Due to the numerous junctions and the highly entangled meshes of city street 
and canals, the random walkers loss their ways and it takes long time for them to 
cross a city roaming randomly along the streets. This can be interpreted as an 
effect of a slight negative pressure involving the random walkers into the city. It 
is obvious that the strength of pressure should vary from one district to another 
within a city and can be different for the different modes of diffusion process 
which have no rigorous bind to the city administrative units being localized on 
the certain groups of streets and canals. We have computed the pressure spectra 
-P(A) forcing the flows of random walkers with eigenmodes A into the compact 
city structures (see Figs. 20-22). Generally speaking, the more junctions a city 
has, the stronger is the drag force: despite Venice has more canals than in 
Amsterdam, the number of junctions between them is less than in Amsterdam, 
and therefore the negative pressure in the latter city is stronger (Fig. 20). The 
numbers of streets in Rothenburg and in the downtown of Bielefeld are equal, 
but there are more crossroads in Bielefeld than in Rothenburg (Fig. 21). 

One can see that pressure profiles have maxima close to A = 1 (correspondent 
to a minimal drug force) for Amsterdam, Venice, and Manhattan. It calls for 
the idea of a "transparency corridor", i.e. a sequence of streets and junctions 
(on which the relevant eigenmodes are localized) along which the city can be 
crossed at minimal time. In Manhattan, the pressure profile is almost zero at 
A = 1. The multiply degenerated eigenmode A = 1 appears in Manhattan due 
to numerous junctions of low accessible streets to Broadway and FDR Drive. 
In Venice, the minimal pressure is achieved on the Grand Canal and Giudecca 
Canal, and it is due to Het IJ and Amstel river in Amsterdam. 

In Tab. 2,we have sketched the figures for the internal energy, entropy, the 
free energy, and heat capacity for all studied cities at = 1 (for the usual 
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random walks model), in purpose of comparison. We also gave the location of 
maxima for the heat capacity profiles. 

6 Discussion and Conclusion 

We have studied the finite Markov chain processes defined on the dual graphs 
of compact urban patterns. The traffic of random walkers, indeed, do not de- 
scribe the actual traffic conditions in the city patterns that we have analyzed, 
but concerns their topological properties giving a sense to the notions of accessi- 
bility and modularity. Our approach can be readily used in order to investigate 
connectedness and efficiency of transportation lines in different complex net- 
works. 

The methods developed in graph theory and in probability theory give us 
a detailed picture of local and global properties of city structures. Dynamical 
modularity representing the localization of dynamical modes on certain streets 
and canals cannot be detected neither from the adjacency matrix of a graph nor 
from the transition probability matrix. We have studied the random walks of 
a large number of massless particles introduced in the dual graphs of compact 
cities and investigated their properties by means of statistical mechanics. The 
spectrum of Laplace operator is broadly distributed for the compact city pat- 
terns and has a bell shape, in general. However, it turns into a sharp peak if 
there is a number of low accessible streets branching of a prime street in a city. 

To detect the dynamical modularity precisely, we have computed the linear 
correlation coefficients between the lists of eigenvector components for all pairs 
of streets (canals) in the cities. Tuning the sensitivity threshold of correlations, 
one can detect the various modules of essentially correlated streets. The com- 
plexity of dynamical modularity can be measured by the information quantity. 
The information versus correlation profile is an individual city fingerprint. 

The statistical mechanics description of random walks in the compact cities 
provides us with a rigorous definition of the unique measure and all statistical 
moments which have a definite relation to the entire topological properties of 
complex networks. Due to a complicated topology of city streets and canals, 
the flows of random walkers acquire nontrivial thermodynamical characteristics 
which can be implemented in studies of interactions between different city mod- 
ules and of city stability with respect to occasional accessibility problems that 
would be of vital importance for the city traffic conditions in urgency. 

The obvious advantage of the method is that the information on the role of 
a given street or canal in the entire city network can be interpreted by the very 
end, and all city modules can be named street by street. 

The method can be generalized for the complex networks which can be de- 
scribed by directed or weighted graphs. For instance, one can consider a graph 
of city plan in which the weights of edges are the actual lengthes of streets. 
In the case of weighted adjacency matrix, the spectrum of the relevant Laplace 
operator acquires the pairs of complex conjugated eigenvalues. Then the compo- 
nents of eigenvectors are also complex, and the coefficients of linear correlations 
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should be computed for the real and imaginary parts of spectrum separately. 
Other approach could be used while studying directed graphs for which the ad- 
jacency matrix is not symmetric. For the random walks defined on the directed 
graphs, the probability that a random walker enters a node is not equal to the 
probability it leaves the node. The ensembles of random walkers introduced on 
the directed graphs can be described by Nelson stochastic mechanics [51], [52] , 
The stationary configurations of random walkers in such models can be studied 
by means of biorthogonal decomposition 53J . 
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Table 1: Some features of dual graphs of the studied city patterns 





\v\ 


\E\ 


diam(G) 


T 


Rothenburg ob.d.T. 


50 


115 


13 


545 


Bielefeld downtown 


50 


142 


14 


551 


Amsterdam canals 


57 


200 


11 


849 


Venice canals 


96 


196 


14 


1550 


Manhattan 


355 


3543 


17 


4557 



Caption for Table 1: Some features of dual graphs of compact city pat- 
terns: \V\ is the number of streets, \E\ the number of junctions (crossroads), 
the graph diameter diam(G) is the maximal graph-theoretical distance between 
any two vertices of the dual graph. In the last column, the the random target 
access time r, the expected number of steps required to hit a node randomly 
chosen in the city from the stationary distribution w. The random target access 
time is independent of the starting point. 
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Table 2: Thermodynamical properties of ensembles of random 
walkers in the compact cities at /3 = 1 (the usual random walks). 





E(J3 = 1) 


509=1) 


F(J3 = 1) 


c 

^max 




Rothcnburg ob.d.T. 


0.946604 


3.885109 


2.93850 


1.60468 


11.1 


Bielefeld downtown 


0.957324 


3.89042 


2.93309 


2.04195 


9.90 


Venice canals 


0.95082 


4.53962 


3.58880 


1.92003 


27.4 


Amsterdam canals 


0.898409 


3.88836 


2.98996 


1.3353; 1.1094 


8.6; 48.75 


Manhattan 


0.990083 


5.86713 


4.87705 


5.87368 


13.65 



Caption for Table 2: E is the internal energy, S is the entropy, F is the 
free energy, C max is the maximal values of heat capacity, and /3c max is the points 
they are achieved. 
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Figure 1: The logarithm of fractions of streets lnn(s) via the logarithm of 
number of their junctions Ins in Manhattan (shown by points). The solid line 
represents for the logarithm of the relevant cumulative distribution In P c {s) [7]. 
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Figure 2: The distributions of mean access times h to the streets in Manhattan 
(dashed line) and Rothenburg ob der Tauber (solid line). The distribution ah 
for Rothenburg o.d.T. exhibits a local maximum at relatively long access times 
( 300 steps) indicating the presence of a number of low accessible streets in the 
town. 
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Figure 3: The distributions of mean access times h to the canals in Amsterdam 
(dashed line) and Venice (solid line). 
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b. 



Figure 4: The city map of Bielefeld downtown (a) is presented together with 
its 3D representations of the dual graph (b). The "A" -part keeps its original 
structure (founded in XIII- XIV cs.); the part "B" which had been redeveloped 
in the XlX-th century. The (xi, yi, Zi) coordinates of the i— th vertex of the dual 
graph in three dimensional space are given by the relevant i— th components of 
three eigenvectors , ,and of the adjacency matrix A of the dual graph 




Figure 5: The distributions of mean access times h to the streets located in 
the medieval part " A" starting from those located in the same part of Bielefeld 
downtown, from "A" to "A" (solid line). The dashed line presents the distri- 
bution of mean access times to the street located in the modernized part "B" 
starting from the medieval part "A" (from "A" to "B"). In average, in takes 
longer time to reach the streets located in " B" starting from " A" . 
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Figure 6: The distributions of mean access times h to the streets located in the 
" B" part starting from " B" (solid line) . The dashed line presents the distribution 
of mean access times to the street located in the " A" part starting from " B" . 
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Figure 7: The eigenvalues of the Laplacian ([3]) defined on the dual graph of 
Bielefeld together with the continuous approximation of its density. 
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Figure 8: The eigenvalues of the Laplacian ([3]) defined on the dual graph of 
Amsterdam canal network together with the continuous approximation of its 
density. 
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Figure 9: The eigenvalues of the Laplacian ([3]) defined on the dual graph of 
Manhattan together with the continuous approximation of its density. 
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Figure 10: a. Plan of the Chelsea village in Manhattan and the matrix plot of 
its nodal domains (b) . All eigenvectors localized on the nodal domains displayed 
in white (black) have always the positive (negative) sign. 
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Figurc 11: Quantity of information (bits) encoded by the dynamical modu- 
larities of Rothenburg and Bielefeld via the threshold of essential correlations, 

< K < 1. 



33 



o.i2; 



i "bits" 



0.1 



0.06 



0.04 



-I \'r- 



0.02 



Venice 




Amsterdam 



K 



Figure 12: Quantity of information (bits) encoded by the dynamical modulari- 
ties of the Venice and Amsterdam canal networks via the threshold of essential 
correlations, < k < 1. 
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Figure 13: Quantity of information (bits) encoded by the dynamical modulari- 
ties of Manhattan via the threshold of essential correlations < k < 1 . 
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Figure 14: The grows of "internal energy" in the models of lazy random walks 
with "temperature" /3 _1 (the inverse parameter of lazy random walks) for the 
canal networks of Amsterdam and Venice. 
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Figure 16: The grows of "internal energy" in the models of lazy random walks 
with "temperature" /3~ 1 (the inverse parameter of lazy random walks) for Man- 
hattan and the theoretical example of a "village" extended along the only prin- 
cipal street. 
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Figure 17: The entropy curves via the inverse parameter of lazy random 
/3 _1 , for the canal networks of Venice and Amsterdam 
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Figure 20: The comparison of pressure spectra P(X) acting on the flows 
random walkers with eigenmodes A into Amsterdam and Venice. 
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Figure 21: The comparison of pressure spectra P(X) acting on the flows of 
random walkers with eigenmodes A into Bielefeld and Rothenburg. 
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Figure 22: The comparison of pressure spectra P(A) acting on the flows of 
random walkers with eigenmodes A into Manhattan and Venice. 
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